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Abstract 



(^ We simulate the assembly dynamics of icosahedral capsids from 

^ subunits that interconvert between different conformations (or quasi- 

.^ equivalent states). The simulations identify mechanisms by which sub- 

C<^ units form empty capsids with only one morphology, but adaptively 

'^ assemble into different icosahedral morphologies around nanoparticle 

t^^ cargoes with varying sizes, as seen in recent experiments with brome 

rt: mosaic virus (BMV) capsid proteins. Adaptive cargo encapsidation re- 

(^ quires moderate cargo-subunit interaction strengths; stronger interac- 

^ tions frustrate assembly by stabilizing intermediates with incommen- 

surate curvature. We compare simulation results to experiments with 
cowpea chlorotic mottle virus empty capsids and BMV capsids assem- 
C^ bled on functionalized nanoparticles, and suggest new cargo encapsi- 

dation experiments. Finally, we find that both empty and templated 
capsids maintain the precise spatial ordering of subunit conformations 
seen in the crystal structure even if interactions that preserve this ar- 
rangement are favored by as little as the thermal energy, consistent 
with experimental observations that different subunit conformations 
are highly similar. 



*To whom correspondence should be addressed: hagan@brandeis.edu 



Introduction 

During the replication of many viruses, hundreds to thousands of protein sub- 
units assemble around the viral nucleic acid to form a protein shell called a 
capsid. In vitro studies show that capsid proteins can form particular empty 
capsid structures with high fidelity; ^ yet capsids adopt different morphologies 
when challenged with nucleic acids ^ or other cargoes ^'^'^ with sizes that are 
not commensurate with the preferred capsid structure. No proposed dynam- 
ical mechanism simultaneously explains precise assembly of empty capsids 
and adaptable encapsidation of cargoes. Understanding how viral compo- 
nents selectively assemble into the structure required for infectivity could 
spur the development of antiviral therapies that block or alter assembly. At 
the same time, engineered structures in which viral capsids assemble around 
synthetic cargoes show great promise as delivery vehicles with adaptable sizes 
for drugs or imaging agents,^ and as subunits or templates for the synthesis 
of nanomaterials with exquisitely controlled sizes and morphologies.^ Real- 
izing these goals, however, requires understanding how properties of cargoes 
and capsid proteins dynamically determine the size and morphology of an 
assembled structure to enable adaptable assembly. 

In this work, we explore the interplay between cargo size and the mor- 
phology of icosahedral capsids with coarse-grained models that describe both 
the dynamic encapsidation of functionalized nanoparticles and the assembly 
of empty capsids. Through our simulations, we uncover a mechanism by 
which subunits faithfully assemble into empty capsids with a single icosahe- 
dral morphology, but also reproducibly assemble into different morphologies 
around nanoparticles with varying sizes, as seen in recent experiments. ^ The 
model predicts that adaptability to cargo size is nonmonotonic with respect 
to the strength of subunit-cargo interactions. This prediction can be tested 
in nanoparticle-capsid assembly experiments by varying the functionalized 
surface charge density on nanoparticles.® 



Assembly of icosahedral viruses. While at most 60 identical subunits can 
be arranged with icosahedral symmetry, Caspar and Klug showed that mul- 
tiples of 60 proteins can form icosahedral capsids, if individual proteins take 
slightly different, or quasi-equivalent, conformations.^ These quasiequivalent 
conformations break the local 3-fold symmetry of the icosahedral face but, 
by assembling with precise spatial ordering of conformations, preserve the 
global icosahedral symmetry of the capsid. Despite their different geometry. 



the proteins interact with each other by interfaces that are substantially sim- 
ilar across different conformations. A complete capsid is comprised of GOT 
subunits, where T is the number of distinct protein conformations (see 2). 

Although recent experiments ^° have begun to characterize sub unit con- 
formations during assembly and equilibrium theories have led to important 
extensions of quasi-equivalence, ^^ the process by which the appropriate quasi- 
equivalent conformations are chosen during assembly remains poorly under- 
stood. Berger and coworkers ^^ showed that Caspar-Klug structures result if 
assembly follows "local rules", in which only subunits with the conforma- 
tion dictated by adjacent subunits can bind to an assembling capsid. There 
are two experimental observations, however, that seem difficult to rational- 
ize with conformation-dependent interactions. (1) How can subunit-subunit 
binding be conformation-specific for viruses in which subunit binding inter- 
faces show little variation between conformations in capsid crystal structures 
(e.g. see 3 and Ref. ^^)? (2) How can subunits that assemble with conforma- 
tional specificity adapt to form capsids with different icosahedral morpholo- 
gies around commensurate cargoes? For example, Dragnea and coworkers ^'^^ 
have demonstrated that brome mosaic virus (BMV) proteins assemble into 
T=l, pseudo-T2, and T=3 capsids around functionalized nanoparticle cores 
with different diameters that are functionalized with carboxylated polyethy- 
lene glycol. 

We explore both of these questions here with a model for the assembly of 
T=l and T=3 capsids from subunits that can interconvert between different 
conformations, with which we simulate the spontaneous assembly of empty 
capsids and the encapsidation of nanoparticles that template assembly of 
T=l or T=3 capsids. By systematically varying the extent to which binding 
between subunits depends on conformation, we show that even a weak con- 
formational dependence (~ ksT) enables robust assembly. In addition, we 
find that an intrinsic bias for subunits to adopt particular conformations, as 
suggested by recent experiments, ^^ can control which icosahedral morphology 
is favored. We show that requiring the model to reproduce experimental ob- 
servations for both empty and full capsids places tight constraints on model 
parameters, and thereby enables insights about morphology control in both 
systems. In particular, we find a narrow, but physically reasonable, range 
of parameters for which only T=3 empty capsids assemble, but T=l capsids 
form on a commensurate nanoparticle. 



Model 

We extend a class of models for T=l capsids^^'^^ (see Refs. ^''' for related 
models), in which subunits have spherically symmetric excluded volume 
and short-ranged, directional attractions between complementary interfaces. 
These interfaces are represented as 'bond-vectors' that rotate rigidly within 
the subunit frame of reference; we design their geometry such that the sub- 
units tile an icosahedron according to capsid crystal structures. Thus, the 
lowest energy collective configurations correspond to "capsids" with GOT 
monomers in a shell with icosahedral symmetry. 

We focus on models for T=l and T=3 capsids, for which the bond vector 
geometries are based on crystal structures of T=l^^ and T=3^^ BMV cap- 
sids. Each subunit represents a protein dimer, the basic assembly unit for 
BMV. ^° The relative positions and conformations of subunits in a capsid are 
determined by associating each two-fold or quasi-two-fold dimer interface 
with a subunit center, as depicted in 2. The orientations of bond-vectors 
and their complementarity are then determined from the relative locations 
of neighboring dimers; as shown in 2 each interface between neighboring sub- 
units is associated with a pair of complementary bond vectors. The resulting 
internal coordinates and list of complementary bond vectors are specified in 
the SI. 

Pair interaction. The attractive interaction between two complementary 
bond-vectors b, and bj (see 1) is minimized when (1) the distance between 
the bond- vectors r^- is minimized, (2) the angle between them, 6^-, is mini- 
mized and (3) the dihedral angle, 0^-, calculated from two secondary bond- 
vectors which are not involved in the primary interaction (see the SI and^^), 
is minimized. Requirement (3) creates an interaction that resists torsion and 
therefore enforces angular specificity in the plane perpendicular to the bond 
vector. The potential is given by equations (1) through (5) 



(a) (b) 



^ h?V--^^ ^^ ^fj = arccos(-b, ■ b 



arccos(— bj ■ bj/|bi||bj|) 



Figure 1: Definition of parameters used to calculate the interaction energy 



U 



between complementary bond vectors (2) and (3). The dihedral angle 
giving the rotation around the plane perpendicular to Rij is not shown. The 
bond vectors here are shown projected onto a plane even though the subunit 
geometries have intrinsic curvature (see 2). 



U = Ure, + Y.U,u{r\^)S{e\^,ct>\^) (1) 



b 
Urep = e(i?-2g) {C{R) + l) (2) 

Uatt = e(r,,-r,Kx-,(>C(r., + 2i)-/:(re)) (3) 

(cos(7r0/0e) + 1) (cos(7r0/0,) + 1) (4) 

C{x) = 4(x-^'-x-^) (5) 

where the index b sums over pairs of complementary bond vectors, 0(x) is 
the Heaviside step function, R is the subunit center-to-center distance, C{x) 
is a 'Lennard- Jones function', and the cutoff values are Tc = 2.5, ^c = 1 and 
0c = TT throughout this work. Throughout this work, lengths have units of 
a, the subunit diameter, energies have units of k^T and times have units of 
^0 = ^^/-D, where D is the subunit diffusion constant. Concentrations are 
defined as Cq = Na^/L^ with A^ subunits and box side length L. 

Conformation dependence of subunit- subunit binding energies. We follow 
convention^ by labeling the different protein (monomer) conformations found 
in the BMV or CCMV crystal structure as A, B, and C A T=3 capsid is 



comprised of 30 CC and 60 AB dimer subunits, while 30 AA dimer subunits 
comprise a T=l capsid; the structure of an A protein monomer visible in the 
T=3 CCMV capsid crystal structure is virtually identical to the structure 
of monomers in the T=l crystal structure. ^^ For simplicity, we assume that 
dimer configurations not present in the crystal structures (e.g. AC, BB) are 
highly unfavorable and thus do not occur. 

If the subunit binding interactions have strict conformational specificity, 
then subunits i and j can only bind via a particular interface b ONLY if the i 
and j have conformational states that interact via that interface in the crystal 
structure (see 2). As noted above, however, there is not strong evidence to 
support such strict conformational specificity. Therefore, we take xtj = 1 
(Eq. 3) if the conformations of subunits i and j are found in the crystal 
structure for interface b, and x'ij = Xs otherwise, where the promiscuity 
parameter Xs varies from for strict conformational specificity to 1 for no 
conformational specificity at all. 

To capture the ability of BMV to form both T=l and T=3 structures, we 
allow subunits to change conformation during the simulation. For simplic- 
ity, we model the transitions as discrete events with no intermediate states, 
implying that subunit conformations are separated by an activation barrier. 
These moves are accepted according to the Metropolis condition: 

P = min{l,exp(-Af/- Ae)} (6) 

where AU is the change in the interaction energy and Ae is the intrinsic free 
energy difference between conformational states, which might correspond 
to the energies associated with different hinge angles described in Ref. ^^ 
(see 3). We set the energy of AB and the CC energy to log 2 (because 
it is symmetric). We vary the remaining energy, ^aAj between and 2.5. 
For simplicity, we do not consider pseudoT=2 capsids in this work, which 
unlike T=l and T=3 capsids, involve binding interfaces that are not seen in 
infectious CCMV viruses. ^^ 

Core- controlled assembly. 

Motivated by recent experiments in which BMV capsid proteins encap- 
sidate functionalized gold nanoparticles,^^ we introduce a rigid sphere with 
radius Rg at the center of the simulation. The sphere interacts with the sub- 
units via a spherically symmetric Lennard- Jones potential, shifted so that a 
subunit at the surface of the sphere has minimum energy 

Us{r) = eMrcS - r^) {C{r,s) - C{r,)) (7) 



(a)Tl 



(b) T3 





Figure 2: (a) Images of brome mosaic virus (BMV) capsid crystal struc- 
tures (T=l,^^ T=3 /^ )overlaid with icosahedral cages. Subunits with quasi- 
equivalent conformations A, B, and C are blue, red, and green, respectively, 
(b) The model capsid geometries. Each subunit represents a protein dimer, 
subunit sizes are reduced to aid visibility, and colors indicate the conforma- 
tions of proteins within a dimer subunit: silver, AA; red, AB; blue, CC. 



where Ves = r~Rs with r the nanoparticle-subunit center-to-center distance, 
and Ss specifies the strength of the subunit-sphere interaction; we consider the 
range Es/k^T G [6, 12]. We consider nanoparticles with Rg = 1.7 and Rg = 
3.2, which are commensurate with T=l and T=3 model capsids, respectively. 

Dynamics simulations. We evolve particle positions and orientations with 
over-damped Brownian dynamics using a second order predictor-corrector 
algorithm. ^^ We intersperse conformational Monte-Carlo moves with dynam- 
ics integration steps such that, on average, each particle attempts to change 
conformation with frequency * 7c = (40to)~^- When there is a nanoparti- 
cle present we simulate a bulk solution with concentration cq by performing 
grand canonical Monte Carlo moves in which subunits far from the nanoparti- 
cle are exchanged with a reservoir at fixed chemical potential with frequency 
consistent with the diffusion limited rate.^^ Since only a single nanoparti- 
cle is considered in each simulation, interactions between nanoparticles are 
not considered - finite nanoparticle concentrations are considered in Ref. ^^ 
Empty capsid simulations have A^ = 1000 subunits in a box with side length 
L = 22.5a. 



* Assembly behavior appears to be largely independent of 7^, except at extreme values: 
assembly is not adaptable at 7c = 0, and assembly is not productive at 7c = 00 



(a) 



(b) 





Figure 3: (a) Subunits with different conformations in CCMV capsids have 
nearly identical structures, but different angles of rotation about the "dimer 
hinge" , which is the axis through the center of the dimer interface. ^^ One 
sub unit is shown on the left (dark blue), and subunit positions on the right 
are indicated by the ribbon representations of two /5-strands to show the 
range of rotation about the hinge for AB (red, 38°), CC (light blue, 42°), 
and AA (green, 45°). (b) Geometries and intrinsic free energies (or hinge 
energies) for model subunits in different conformations. The only geometrical 
difference between AA and CC subunits is the "hinge" angle between the left 
and right pairs of bond vectors - the AA dimer has higher curvature (into 
the page) than the CC. Bond vectors are depicted with arrows and the angles 
between them are indicated in degrees - the angles do not add up to 360° 
because the vectors do not lie in a plane. 
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Figure 4: The fraction of subunits in well-formed T=l (■),T=3 (•) or both 
(♦) capsids for varying conformational free energy, Sp^x, in empty capsid 
simulations. 



Results 

Empty capsids. In order to understand which system parameters control 
assembly morphology, we performed dynamical empty capsid simulations for 
varying values of the conformational free energy, e^A, and the conformation 
specificity parameter, Xs, which controls the conformational dependence of 
subunit-subunit binding energies. Prior work^^ has shown that assembly 
yields are nonmonotonic in the parameters that control the driving forces 
for assembly, the binding energy, ^b, and the angular tolerance, 6c. For 
empty capsid simulations in this work, we consider only optimal values of 
these parameters, ^b = 12.0, 6^ = 1.0, for which subunits that are capable of 
forming only a single morphology assemble with high fidelity. ^^ 

The conformational free energy controls polymorphism in empty capsid as- 
sembly. As a measure of morphology control, we monitor the fraction of 
subunits in T=l and T=3 capsids, which are defined as connected clusters 
comprising 30 or 90 subunits, respectively, in which each subunit has 4 bonds. 
The yield of each morphology as a function of the conformational free en- 
ergy, ^AAj is shown in 4 for an observation time of 200, OOOto, at which point 
assembly has roughly saturated. We observe a crossover from high yields of 
T=l capsids for ^aa < to predominantly T=3 capsids for eaa > 1-0, with 
mixed morphologies in the intermediate region. 

Although the transition between T=l and T=3 capsids only requires 
a change in the conformational free energy of approximately the thermal 
energy, ksT, the width of the transition is consistent with kinetic rather 
than thermodynamic control of the dominant morphology. As shown in Fig 
6b of the SI, fitting the fraction of subunits in T=l capsids, /ti to the form 
/ti(£^aa) = 1/(1 + 31og2exp[— nc(eAA ~ log2)]), yields a 'critical' size for 
morphology determination of Nc G [6.5,8.0], while an equilibrated system 
would give a much sharper transition with Uc = 30. 

To better understand this result, we estimated the size of a "critical nu- 
cleus" , or an intermediate which is more likely to grow into a complete cap- 
sid then to dissociate into free subunits. A simulation with 10,000 subunits 
was run until many small assembly intermediates (2-12 subunits) assembled. 
From that configuration, further simulations were integrated with different 
random number seeds, and each initial intermediate was tracked in every 
simulation. We then estimate the "commitor probability" for each interme- 
diate as the fraction of simulations in which it grows to completion before 
dissociating. The average commitor probability as a function of intermediate 
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Figure 5: The fraction of subunits in well-formed capsids as a function of the 
conformational specificity parameter Xs- 



size (figure 6a in the SI) suggests that the critical nucleus is 7, which is con- 
sistent with the critical morphology size estimated above. We note, however, 
that this value is only a rough estimate, since the identity of critical nuclei 
depends on additional parameters, such as the number of bonds and closed 
polygons - i.e. the intermediate size alone is not sufficient for a good reaction 
coordinate. 

Faithful assembly requires only weak conformational specificity. As shown 
in 4, varying the intrinsic conformational free energy, ^aa, leads to differ- 
ent morphologies, but does not significantly affect the yields of well-formed 
icosahedrons. As shown in 5, assembly yields are also insensitive to the con- 
formation specificity for Xs ^ 75%, while higher values (indicating less speci- 
ficity) lead to predominantly malformed capsids. These malformed capsids 
are primarily closed but strained structures with disordered arrangements of 
pentamers and hexamers that do not have icosahedral symmetry (see Figs. 
2a and 2b in the SI). 

Core-controlled assembly. To understand the effect of cargo properties 
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on morphology, we simulated assembly in the presence of a model nanoparti- 
cle with varying conformational free energy, ^aa, and core-sub unit attraction 
strength, 5s- Optimal conditions for the nanoparticle encapsidation experi- 
ments have higher pH and thus weaker subunit-subunit binding interactions 
than empty capsid experiments. * Therefore, when simulating assembly with 
a model nanoparticle, we consider subunit-subunit binding parameters of 
Eh = 10.0, 6c = 1.0 for which assembly is favorable on the nanoparticle, but 
no spontaneous assembly in bulk solution occurs. ^^ As a measure of assembly 
effectiveness, we monitor the packaging efficiency, /s, which is defined as the 
fraction of independent trajectories in which a well- formed capsid assembles 
on the nanoparticle. Consistent with prior work, ^^ parameter values that lead 
to a single morphology of empty capsids enable efficient encapsidation of a 
commensurate sphere with ^ 100% efficiency. To understand polymorphism, 
we monitor packaging efficiencies around a T=l sized sphere while varying 
the conformational free energy from values that favor T=l empty capsids 
(^AA ^ 0.25) to values favoring T=3 empty capsids (^aa ^ 1)- The results 
are shown in 6a for several core-subunit interaction strengths. At each ^s, 
there is a relatively sharp crossover from high yields to no successful assem- 
bly. Significantly, for optimal values of Eg, there is a range of 1 < £aa ^ 1-5 
for which spontaneous assembly faithfully produces T=3 capsids, but T=l 
capsids form on the nanoparticle with high efficiency. 

Mechanisms of core- controlled polymorphism. In contrast to assembly 
around a commensurate sphere, ^^ packaging efficiencies do not increase mono- 
tonically with core-subunit interaction strength (Es), as shown for several val- 
ues of ^AA in 6b. Assembly yields that are nonmonotonic with the variation 
of an interaction parameter are a hallmark of competition between thermo- 
dynamics and kinetics. At low 5s, the core-subunit interaction strength is 
not large enough to stabilize a T=l capsid nucleus and so T=l assembly is 
either thermodynamically unfavorable or has an insurmountable activation 
barrier. At high Eg, on the other hand, capsid nuclei of any morphology are 
stabile and T=3 partial capsids are common. Beyond a certain size, subunits 
that add to an incommensurate (T=3) partial capsid cannot simultaneously 
interact strongly with the core surface and subunits already in the interme- 
diate (see 7a). For the parameters considered, subunit-subunit interactions 
are too weak to drive significant assembly away from the attractive core sur- 
face, and assembly is frustrated. At high 5s, however, the partial capsid is 
metastable and blocks a significant fraction of the core surface, thus hindering 
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the formation of potentially productive T=l nuclei. 

At optimal strengths of the core-subunit interaction, partial capsids un- 
dergo significant size fluctuations because subunit unbinding is frequent. 
Therefore, many nuclei can form on a given sphere within the observation 
time, until a T=l nucleus grows to completion. 7a shows a metastable 
partial capsid forming on an incommensurate size sphere which eventually 
disassembles, allowing a T=l capsid to grow to completion (7b). 

Time scale for annealing of surface- ads orbed complexes. We simulated 
assembly with a T=l-size nanoparticle over a wide range of subunit concen- 
trations, 3 X 10"^ < Co < 0.06, which corresponds to a range of 7 — lOOO/iM 
for a dimer subunit diameter of 4.2 nm. The analysis in Ref. ^^ suggests that 
nanoparticle systems with different Es and Cq should be compared in terms 
of Csurf) the equilibrium surface concentration of subunits with no assembly^. 
The equilibrium surface concentration can be determined from simulations 
with eb = or by calculating the chemical potential of adsorbed subunits 
(see the SI). The packaging efficiencies for various Eg and cq are plotted as 
a function of Csurf in 9. We see that optimal assembly for all concentrations 
collapses onto the same value of Cgurf ~ 0.4, while higher values of Cgurf niainly 
lead to trapped incommensurate partial capsids. Interestingly, successful as- 
sembly occurs at higher Csurf if the rates of subunit adsorption are decreased 
below the diffusion limited rate by decreasing the frequency of subunit ex- 
changes in the bath (see Figure 6 in the SI), when subunits adsorb more 
slowly compared to the timescale for annealing of surface-adsorbed com- 
plexes. We therefore note that effects which increase the surface-annealing 
timescale, such as barriers to diffusion of strongly adsorbed subunits, could 
further promote frustrated states. 



^Csurf ~ 4^R2 where Ug is the number of subunits with strong (< —2kBT) interactions 
with the nanoparticle. 
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Figure 6: The fraction of trajectories in which a complete capsid assembles 
on a T=l-size nanoparticle as functions of (a) £aa for several Eg and (b) 
^s at several ^aa, All packaging efiic[Lfncies in this work are shown for an 
observation time of 40, OOOto which is 100 times the average assembly time 
under optimal conditions . 
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Figure 7: a) A partial T=3 capsid grows on a T=l-sized nanoparticle; the 
subunits near the sphere are stabihzed by its attraction but the mismatch 
between the capsid and nanoparticle geometry forbids the addition of new 
particles, b) After some time, the metastable partial capsid disassembles 
and a T=l capsid grows to completion. An animation of the entire process 
is available in the SI. 
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Figure 8: Packaging efficiencies as functions of time for encapsidation of T=l- 
size and T=3-size nanoparticles for parameters at which T=3 capsids are the 
lowest free energy empty capsid structures, £:aa = 1-3. Other parameters are 
Es =8.0 and cq = 0.01. 
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Figure 9: Packaging efficiencies on a T=l-size nanoparticle for several bulk 
subunit concentrations, Cq are shown for varying subunit-surface energies, 
plotted as functions of the surface concentration at which subunit adsorption 
saturates, Cgurf- The conformational free energy is ^aa = 1-3. 
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Larger Capsid Morphologies 

To further explore the requirements for adaptabihty to template size, we 
extended our model to include T=4 capsids, which requires at least one new 
subunit geometry in addition to AB and CC (see page 3 of the SI). We note, 
however, that Sun et al.^ observe only disordered structures on T=4-size 
nanoparticles, which may suggest T=4 subunit geometries are not favorable 
for BMV proteins. As shown in figure 4 in the SI, these model subunits 
are capable of forming T=3 or T=4 capsids around nanoparticles ranging 
from T=3 to T=4 size; however, yields of perfect capsids are generally lower 
( rS 50%) and more sensitive to parameter values than for the T=l/T=3 case. 
In addition to intermediates that cannot close around an incommensurate 
core (as described above), there are significant yields of closed but disordered 
asymmetric shells. This result is consistent with the fact that the difference 
in curvature between T=3 and T=4-size nanoparticles is small (~ 12%) and 
thus templating for the commensurate capsid geometry is weak, especially 
for small intermediates. 

Discussion 

Estimating the conformational free energy from experiments. Comparison 
of model predictions to experimental observations of capsid morphology sug- 
gests a potential correspondence between ranges of the conformational free 
energy, ^aA) and certain systems. In particular, CCMV and BMV capsid 
proteins assemble into exclusively T=3 empty capsids under conditions com- 
monly employed in vitro, while deletion of the N-terminal residues from the 
proteins of either virus enables assembly of T=l capsids. ^^ In particular, 
proteins of the A^A34 CCMV mutant, with 34 N-terminal residues deleted, 
assemble into a mix of T=3 (~ 10% of assembled material), T=l (^ 40%) 
and heterogeneous pseudoT2 assemblages (~ 50%), could correspond to the 
range 0.3 ^ £aa ~ 0.8 for which polymorphism is observed. We note that this 
range could shift somewhat, however, if there are different binding energies 
for different complementary interfaces and because we have not considered 
pseudoT2 capsids in this work, since they require interfacial contacts that 
are not seen in the T=3 crystal structure (unlike T=l capsids). ^^ 

The observation that wild-type CCMV and BMV proteins form exclu- 
sively T=3 capsids can only suggest a lower bound for the wild-type confor- 
mational free energy, ^aa ~ 1-2, but additionally considering nanoparticle 
experiments allows a more precise estimate. Comparison of Figs. 4 and 6b 
identifies only a narrow range 1.1 ^ ^aa ~ 1.5 for which predominantly 



T=3 empty capsids form, but T=l capsids efficiently encapsidate T=l-sized 
nanoparticles, as observed in BMV-nanoparticle experiments.^'* Although 
CCMV proteins have not been used in nanoparticle experiments, pseudo-T2, 
T=3, and larger (but asymmetric) CCMV capsids assemble around inorganic 
polyelectrolytes^ and nanoemulsion droplets.* 

Although we have focused on the relationship between nanoparticle sur- 
face charge and assembly effectiveness, other model parameters can be var- 
ied in experiments as well. For instance, decreasing pH (or increasing ionic 
strength) can increase sub unit binding free energies^* (sb) while simultane- 
ously decreasing the subunit-nanoparticle interaction strength (ss), since the 
nanoparticle surfaces are functionalized with carboxylated PEG. Thus, the 
optimal pH for encapsulation of nanoparticles is larger than the optimal pH 
for empty capsid assembly.^ Therefore we have simulated assembly in the 
presence of nanoparticles with a subunit binding energy of Sb = 10, which 
is lower than the optimal subunit binding energy ^b = 12 for empty capsid 
assembly; the interdependence of assembly effectiveness on subunit-subunit 
and the subunit-nanoparticle interaction strengths is explored in Ref^^. We 
also note that the conformation free energy (hinge energy) could depend on 
pH. 

Implications for quasi- equivalence. The predictions of our model many 
shed light on the mechanisms by which subunits can assemble with precise 
spatial ordering of different conformations even though the bonding inter- 
faces in different conformations are structurally similar. In particular, the 
model predicts that assembly products, and to some extent assembly times, 
are insensitive to the inter-subunit conformational binding specificity, Xs, for 
Xs ~ 25%. For the parameters used in this work, the free energy per bond in a 
complete capsid is approximately 3.5 ksT (this estimate includes the entropy 
penalty for the subunit binding, see*^). Thus successful assembly requires 
only that conformational pairings not seen in the crystal structure differ by 
~ ksT from native pairings, which could arise from only minor variations 
in binding interfaces. Because optimal assembly occurs for weak subunit- 
subunit interactions, when subunit binding is only slightly more favorable 
than unbinding, ^^ a small difference in subunit binding free energies, com- 
bined with the strain caused by the geometrical incompatibilities that result 
from non-native bonding, strongly favors a Caspar-Klug capsid structure. 

Suggested experiments. A key prediction of our work is that assembly on a 
nanoparticle with a size that does not match the lowest free energy capsid 
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morphology is impeded by partial capsids whose curvature is inconsistent 
with the nanoparticle surface. These frustrated states are revealed in sev- 
eral ways that may be accessible to experiments. First, simulated packaging 
efficiencies on T=l-size nanoparticles, for sub units that form T=3 empty cap- 
sids, are nonmonotonic with respect to the nanoparticle-subunit interaction 
strength (6b). This parameter is controlled in experiments by functionalizing 
nanoparticle surfaces with different ratios of anionic and neutral molecules. ^ 
In Ref. ^^ we show that the chemical potential of adsorbed subunits, and hence 
the equilibrium surface concentration Csurf (see 9), can be estimated from the 
surface density of functionalized charge; Cgurf can be varied over the range 
< Csurf ^ 0.8. Although the degree of nonmonotonicity depends on fac- 
tors such as subunit adsorption rates and surface annealing rates, frustration 
can still be observed by differences in assembly kinetics on T=l and T=3- 
size nanoparticles even for conditions in which eventual packaging efficiencies 
reach 100% as in 8. 

Sumniary. We have performed simulations with a model for assembly 
of subunits into empty capsids and around nanoparticles that template the 
assembly of different morphologies. The simulations uncover mechanisms by 
which assembly can adapt to form different morphologies when challenged 
with a template that does not match the preferred empty capsid structure. 
Predicted assembly pathways include frustrated partial capsid intermediates 
with curvatures that do not match the template, which leads to predicted 
differences in assembly kinetics and effectiveness on nanoparticles with dif- 
ferent sizes, that can be tested in experiments. These findings may shed light 
on the role of nucleic acids in assembly during viral replication, and demon- 
strate that the interplay between the geometries of different components is 
an important consideration for the design of nanostructured materials. 

Outlook. Extensions to this study could include an explicit represen- 
tation of intra-subunit degrees of freedom (i.e. hinge motions of the pro- 
tein dimer), and template degrees of freedom. As suggested by an anony- 
mous reviewer, a flexible spherical template could represent a nucleic acid 
molecule with signiflcant structure due to base pairing. We flnd that assem- 
bly around a flexible spherical template has qualitatively similar results to 
those reported here (for some ranges of sphere flexibility), although the opti- 
mal sub unit-template interaction strength increases somewhat; these results 
will be presented in a future publication. 
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In order to focus on the effects of template-capsid geometry mismatches 
on assembly, we have not considered other potential sources of frustration, 
such as impeded diffusion for subunits that interact strongly with the nanopar- 
ticle surface. The coupling of multiple sources of frustration could have in- 
teresting consequences. 

Supplementary Materials 



Supplementary materials are available online. 
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